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Abstract 

Because of Pauli suppression effects the Nir decay mode of the free 
A is not of importance in hypernuclei with A > 10. Rather the decay 
of such hypernuclei proceeds via the nucleon-stimulated mode AN — ► 
NN, analysis of which presents a considerable theoretical challenge 
and about which there exists only a limited amount of experimental 
information. Herein we confront existing data with various theoretical 
analyses which have been developed. 



1 — Introduction to Hypernuclear Decay 

The properties of the lambda hyperon are familiar to all of us. Having a 
mass of 1116 MeV, zero isospin and unit negative strangeness, it decays 
nearly 100% of the time via the nonleptonic mode A — ► Ntt and details can 
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be found in the particle data tables U 



rA " 263pi RR - A ^{ n7T o 357% (1) 

The decay can be completely described in terms of an effective Lagrangian 
with two phenomenological parameters 

H w = g w N{l + Kj 5 )T-$ 7r A (2) 

where g w = 2.35 x 1CT 7 , k = —6.7 and A is defined to occupy the lower entry 
of a two component column spinor. The kinematics are such that for decay at 
rest the final state nucleon emerges with energy about 5 MeV, which means 
that the corresponding momentum is p^ = \/2M N E N « 100 MeV. 

Now, however, consider what happens if the Lambda is bound in a hyper- 
nucleus. |2| In this case, even neglecting binding energy effects, the 100 MeV 
momentum of the outgoing nucleon is generally much less than the Fermi 
momentum of the nucleus so that the decay will be Pauli blocked. A very 
simple estimate of this effect can be generated within a simple Fermi gas 
model, wherein, neglecting any effects of binding energy or of wavefunction 
distortion, one finds 



1 r - = 1 -\Y. N nAnjZ\h(Kr)\lS±) (3) 

njl 



r A n 2 



with N n j£ being the occupation number for the indicated state. The result 
of this calculation reveals that the importance of such pionic decays rapidly 
falls as a function of nuclear mass — 

J_ f A=10 A = 25 ... 

r A 1 1/20 1/120 ... K) 

However, while the existence of the nuclear medium suppresses the Ntt mode, 
it also opens up a completely new possibility, that of the nucleon-stimulated 
decay — AA^ — > NN. Assuming that the energy is shared equally between 
the outgoing pair of nucleons one has then En ~ \{m^ — tun) ~ 90MeV. 
The corresponding momentum is pn ~ 400MeV and is well above the Fermi 
momentum, so that Pauli suppression is not relevant. According to the 
above arguments the importance of this nonmesonic (NM) mode compared 
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Figure 1: Calculated ratio of pionic hypernuclear decay to free lambda decay 
rate. 
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to its mesonic counterpart should rapidly increase with A, and this expec- 
tation is fully borne out experimentally, as shown in Figure 1. A theory 
of hypernuclear weak decay then has basically nothing to do with the pio- 
nic mode favored by a free A and must deal with the much more complex 
AN — > NN process. The observables which can be measured experimen- 
tally and should be predicted by theoretical analysis include 

i) the overall decay rate T NM ; 

ii) the ratio of proton-stimulated (Ap — > np) to neutron-stimulated (An — > 

nn) decay— T p NM /T n NM = T NM (p/n); 

iii) the ratio of parity-violating to parity-conserving decay — 

^n^m/^nm = Tnm(PV/PC) 

— which is measured, e.g., via the proton asymmetry in polarized hy- 
pernuclear decay 

iv) final state n,p decay spectra; 

v) etc. 

The present experimental situation is somewhat limited. Most of the early 
experiments in the field employed bubble chamber or emulsion techniques. It 
was therefore relatively straightforward to determine the ratio of the decay 
rates of the two modes, but much more difficult to measure the absolute 
rates. This changed when an early Berkeley measurement on fO yielded the 
value U 

^3 ±1 (5, 

1 A 

However, this was still a very low statistics experiment with sizable back- 
ground contamination. Recently a CMU-BNL-UNM-Houston- Texas- Vassar 
collaboration undertook a series of direct timing — fast counting — hypernuclear 
lifetime measurements yielding the results summarized in table 1.0 In addi- 
tion, there exist a number of older emulsion measurements in light (A < 5) 
hypernuclei, details of which can be found in a recent review article. [0] How- 
ever, the only experimental numbers for heavy systems are obtained from 
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12/~t 
A U 


1 r 

A 


0.41 ±0.14 


1.25 ±0.16* 


1.14 ±0.20 


T NM (p/n) 


1.07 ±0.58 


0.96l£* 


u - ' °-0.35 



Table 1: Experimental BNL data for nonmesonic hyperon decay. *Note that 
we have scaled the experimental number to account to exclude the pionic 
decay component. 



delayed fission measurements on hypernuclei produced in p-nucleus collisons 
and are of limited statistical precision 

T(23 8 U) = (1.0±J;g) x 10- 10 sec. T(f 9 Bi) = (2.5±?:g) x 10- 10 sec. (6) 



The problem of dealing with a weak two-body interaction within the nu- 
cleus has been faced previously in the context of nuclear parity violation, 
and one can build on what has been learned therein. || Specifically, the weak 
interaction at the quark level is shortranged, involving W,Z-exchange. How- 
ever, because of the hard core repulsion the effective NN effects are mod- 
elled in terms of long-range one-meson exchange interaction, just as in the 
case of the conventional strong nucleon-nucleon interaction, |§ but now with 
one vertex being weak and parity-violating while the second is strong and 
parity-conserving. The exchanged mesons are the lightest ones — 7r =l= , p, u — 
associated with the longest range. (Exchange of neutral spinless mesons is 
forbidden by Barton's theorem. |10||) 

A similar picture of hypernuclear decay can then be constructed, but with 
important differences. While the basic meson-exchange diagrams appear as 
before, the weak vertices must now include both parity-conserving and parity- 
violating components, and the list of exchanged mesons must be expanded 
to include both neutral spinless objects (7r°,?7 ) as well as strange mesons 
(K, K*), as first pointed out by Adams. fTT| Thus the problem is considerably 
more challenging than the corresponding and already difficult issue of nuclear 
parity violation. 

One of the significant problems in such a calculation involves the evalua- 
tion of the various weak amplitudes. Indeed, the only weak couplings which 
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Figure 2: Meson exchange picture of nuclear parity violation. 



are completely model-independent are those involving pion emission, which 
are given in Eqn. 2. In view of this, a number of calculations have included 
only this longest range component. Even in this simplified case, however, 
there is considerable model-dependence, as the results are strongly sensitive 
to the short-ranged correlation function assumed for the nucleon-nucleon in- 
teraction, as will be seen. Below we shall review previous theoretical work 
in this area and detail our own program, which involves a systematic quark 
model- (symmetry-) based evaluation of weak mesonic couplings to be used 
in hypernuclear decay calculations. 



2 — Hypernuclear Decay in Nuclear Matter 

As discussed above one of the significant problems in the calculation of hy- 
pernuclear decay involves the evaluation of the various weak NNM vertices. 
Indeed, the only weak couplings which are completely model-independent are 
those involving pion emission, which are given in Eqn. 2. In view of this, 
a number of calculations have included only this longest range component. 
Even here there is considerable model-dependence, however, as the results 
are strongly sensitive to the short-ranged correlation function assumed for 
the nucleon-nucleon interaction. As a warm-up to a realistic calculation then 
we can begin with a pion-exchange-only calculation in "nuclear matter" (i.e a 
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Transition 


Operator 


'So - %(/ = 1) 


|a(g 2 )(l - a A • a N ) 


% - 3 P (7 = 1) 


\b{q 2 )(a A - a N ) ■ q(l - a A ■ a N ) 


^-3^(7 = 0) 


\c{q 2 ){Z + a-a N ) 


3 SW 3 7A(7 = 0) 


^d(q 2 )(a A ■ qa N ■ q - |<7 A • oW) 


3^^i Pl(/ = ) 


^fe(q 2 )(a A - o^) • g(3 + a K ■ a N ) 


3 5l ^3p i(/=1) 





Table 2: Transition operators of allowed AiV — > NN transitions from relative 
S-states. Here q specifies the relative momentum of the outgoing nucleons 
while a a, <Jn operate on the AN, NN vertices respectively. 

simple Fermi gas model with N n = N p and Pf ~ 270 MeV) with and without 
nucleon-nucleon correlation effects. Here the A — A" relative momentum is 
very soft so that only 1 S'o and 3 Si initial states are assumed to be involved. 
Then 

T NM = J d 3 h J d 3 k 2 J J d 3 ki5 4 (pi-p f ) 

x ^£|</|?C|i>| 2 = £W/?<-a) (7) 

spin a/3 

We can break this down further by identifying effective transition oper- 
ators for the various partial wave channels which contribute to the decay 
process — cf. Table 2 — in terms of which we find for the total nonmesonic 
hypernuclear decay rate 

T NM = / p 2 dpqm N (\a\ 2 + \b\ 2 + |c| 2 + \d\ 2 + \e\ 2 + 3|/| 2 ) (8) 

onfi AN Jo ' 

where u,\ m = — ^ — arises from the switch from the nuclear rest frame to 
the AA" center of mass frame and p,q are related by 

p 2 

q 2 = m N (m A - m N ) + (9) 

The results obtained by various groups are displayed in Table 3. 
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Adams Q 


McK-Gib I| 


Oset-Sal 13 


UMass nj 


r ^r A r M (no corr.) 


0.51 


4.13 


4.3 


3.84 


^T NM {C0TT.) 


0.06 


2.31 


2.1 


1.82 



Table 3: Non-mesonic hypernuclear decay rates calculated by various groups 
using pion-exchange only in "nuclear matter." 





Adams] 11 


McK-Gib [ 12 1 


Oset-Sal 13] 


UMass 


^NMip/n) (no corr.) 


19.4 






11.2 


T NM {p/n) (corr.) 


2.8 






16.6 



Table 4: Proton to neutron stimulated decay ratios for pion-only exchange 
in "nuclear matter." 

Obviously there is basic agreement except for the pioneering calculation 
of Adams. The problems with his calculation were two — Adams used an 
incorrect value of the weak coupling constant g w as well as too-strong a tensor 
correlation, both of which tended to reduce the calculated rate. When these 
are corrected the corresponding numbers are found to be 3.5 (no correlations) 
and 1.7 (with correlations) and are in basic agreement with other predictions. 
From this initial calculation then we learn that the basic nonmesonic decay 
rate is indeed anticipated to be of the same order as that for the free A and 
the important role played by correlations. 

A second quantity of interest which emerges from such a calculation is 
the p/n stimulated decay ratio, given by 



T NM (p/n) = 



jWF p 2 dpg (\ a 


2 +|&| 2 + 2 


|c 


2 + 3|d| 2 + 3 


e 


2 + 3|/| 2 ) 


Jo 


p 2 dpq(2 \a\ 


2 + 2 |6|2 + 6|/|2) 



(10) 



and which has been calculated by two of the groups, yielding the results 
shown in Table 4. An interesting feature here is that the numbers come 
out so large — proton stimulated decay is predicted to predominate over its 
neutron stimulated counterpart by nearly an order of magnitude. The reason 
for this is easy to see. In a pion-exchange-only scenario the effective weak 



8 



interaction is of the form 



H w ~ gNfN ■ NtA (11) 

Then An — > nn ~ g but Ap — > np ~ (— 1 — (\/2) 2 )g = —2>g since both charged 
and neutral pion exchange are involved. In this naive picture then we have 
^NM{p/n) ~ 9, in rough agreement with the numbers given in Table 4. 

Armed now with theoretical expectations, we can ask what does experi- 
ment say? The only reasonably precise results obtained for nuclei with A > 4 
are those measured at BNL on |He, ^ 2 C and j^B, which are summarized 
in Table 1. We observe that the measured nonmesonic decay rate is about a 
factor of two lower than that predicted in Table 3 while the p/n stimulation 
ratio differs by at least an order of magnitude from that given in Table 4. 
The problem may be, of course, associated with the difference between the 
nuclear matter within which the calculations were performed and the finite 
nuclear systems which were examined experimentally. Or it could be due to 
the omission of the many shorter range exchanged mesons in the theoretical 
estimate. (Or both!) 

Before undertaking the difficult problem of finite nuclear calculations, it 
is useful to first examine the inclusion of additional exchanged mesons in our 
calculations. As mentioned above, a primary difficulty in this approach is 
that none of the required weak couplings can be measured experimentally. 
Thus the use of some sort of model is required, and the significance of any 
theoretical predictions will be no better than the validity of the model. One 
early attempt by McKellar and Gibson, for example, included only the rho 
and evaluated the rho couplings using both SU(6) symmetry methods as well 
as the well known but flawed factorization approach. |12] Well aware of the 



limitations of this method, they allowed an arbitrary phase between the rho 
and pi amplitudes and they renormalized the factorization calculation by 
a factor of 1/ sin 9 C cos 9 C in order to account for the AI = | enhancement . 
Obviously this is only a rough estimate then and this is only for the rho meson 
exchange contribution! A similar approach was attempted by Nardulli, who 
calculated the parity conserving rho amplitude in a simple pole model and 
the parity violating piece in a simple quark picture. |15| Results of these 
calculations are shown in Table 5 

To my knowledge, the only comprehensive calculation which has been 
undertaken to date is that of our group at UMass. In the case of the parity 
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McK-GibH n + p 


McK-Gib[12| tt - p 


Nard.[I| 




1 r 


3.52 


0.72 


0.7 



Table 5: Nonmesonic decay rates in nuclear matter in pi plus rho exchange 
models 



violating interaction we employed a variant of the (broken) SU(Q) W symme- 
try calculations which were employed successfully by Desplanques, Donoghue 
and Holstein to calculate the various weak NNM couplings in the case of nu- 
clear parity violation. In this approach there exist three (in principle un- 
known) reduced matrix elements which, when multiplied by the relevant gen- 
eralized "Clebsch-Gordan" coefficients, relate all such parity-violating am- 
plitudes. Two of these are determined empirically in terms of experimental 
hyperon decay data, while the third is given by a factorization calculation. 
While the success of this approach in the case of nuclear parity violation is not 
without question, [T^] this procedure provides a plausible and unambiguous 
approach to the problem. 

More difficult is the determination of the parity-conserving weak cou- 
plings. In this case we employ a pole model using the diagrams shown in 
Figure 3. What is needed here are the weak parity-conserving amplitudes for 
A — N and £ — N transitions, which we determine via the current algebra 
(chiral symmetry) relations 

lim(7r°n|74-)|A> = ^<n|[i$, «H]|A) = -l-<n|wH|A) 

lim< 7 r°p|fti- ) |£+> = ^ V \[Ft«M w \^) = ^m^) (12) 

and the weak K — n coupling which is similarly given in terms of the exper- 
imental Kirn decay amplitude 

Artt = -iF w ^2 ^7( 7r ° 7r °l^|-^' )phy S ical (13) 

m K — m~ 



Again this procedure has well documented flaws. [17] However, in the present 



context it is reasonable successful and for a first generation calculation, we 
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Figure 3: Pole diagrams used in evaluation of weak parity-conserving AiV 
NN couplings. 



consider it to provide a reasonable estimate for the parity conserving weak 
couplings. 

Combining with the various strong meson couplings we can now substitute 
into the diagrams shown in Figure 3 to generate the many effective parity 
conserving two-body operators which can be used to evaluate the nonmesonic 
decay amplitudes. Details of these procedures are given in ref. 10. Using the 
resultant two-body operators one can then generate the various predictions 
for nonmesonic decay in nuclear matter. Results are summarized in table 6, 
where we specifically identify the contributions from various channels. 

The results are very intriguing. The overall decay rate is reduced some- 
what from its pion-exchange-only value, in agreement with the experimental 
results. More striking is the modification of the p/n ratio and in the ratio of 
parity violating to parity conserving decay, defined as 



values of which are shown in Table 7. We observe that inclusion of additonal 
exchanges plays a major role in reducing the p/n ratio from its pion-only- 





(14) 



f£ AN p 2 dpq(\a\ 2 + |c| 2 + 3|rf| 2 ) 
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7T 


7T 








(no corr.) 


(corr.) 


7T + p 


u, K, K* 


1 c . lc 
oo ~~ *■ OO 


.010 


.000 


.001 


.001 


% - 3 P 


.156 


.037 


.052 


.018 




.312 


.117 


.113 


.456 


3 Si - ^ 


.468 


.128 


.100 


.110 


3 c - 3c 
Ol — ' Ol 


.010 


.789 


.589 


.202 


3 o\ -> 3 P1 


2.93 


.751 


.693 


.444 


Total 


3.89 


1.82 


1.55 


1.23 



Table 6: Decay rates for various combinations of meson exchange in nuclear 
matter. 





T NM (PV/PC) 


TtvmGoM) 


7T (no corr.) 


0.14 


11.2 


7r (with corr.) 


0.18 


16.6 


7T + p 


0.21 


13.1 


n,p,uj,r],K, K* 


0.90 


2.9 



Table 7: The parity violating to parity conserving and p to n ratios for 
hypernuclear decay in "nuclear matter." 
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exchange value. The resulting value of 2.9 is still somewhat larger than the 
experimental values shown in Table 1 but clearly indicate the presence of 
non-pion exhange components. 

The reason that kaon exchange in particular can play such a major role 
can be seen from a simple argument due to Gibson |L8j who pointed out that 
since the final NN system can have either 1=0 or 1=1, the effective kaon 
exchange interaction can be written as 

C e ff = A (pp + nn)nk + Ai(2nppA — (pp — nn)nA) 

~ (A - 3A 1 )ppnA + (A + A x )nnnA (15) 

where the second line is obtained via a Fierz transformation. Since for parity 
violating kaon exchange we have A ~ 6Ai we find||18[| 



Tnm(p/u) = (^^) -1/5 (16) 

which clearly indicates the importance of inclusion of non-pion-exchange 
components in predicting the p/n ratio. 

A second strong indication of the presence of non-pion-exchange can be 
seen from Table 7 in that the rate of parity violating to parity conserving 
transitions is substantially enhanced by the inclusion of kaon and vector me- 
son exchange as compared to the simple pion-exchange-only calculation. We 
can further quantify this effect by calculating explicitly the angular distri- 
bution of the emitted proton in the Ap — > np transition (there can be no 
asymmetry for the corresponding An — > nn case due to the identity of the 
final state neutrons), yielding 

W p {9) ~ l + aP A cos0 (17) 

where 

tf ANkF p 2 dpqf Ref*(V2c + d) 

f£ ANkF pidpql(\a\ 2 + |6| 2 + 3|c| 2 + 3|c/| 2 + 3|e| 2 + 3|/| 2 1 ' 

is the asymmetry parameter. Results of a numerical evaluation are shown 
in Table 8 so that again inclusion of non-pion-exchange components has a 
significant effect, increasing the expected Ap — > np asymmetry by more than 
a factor of two. This prediction of a substantial asymmetry is consistent with 
preliminary results obtained for p-shell nuclei at KEK.flT9f 
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7r-no corr. 


7T-COIT. 


all exch. 


a 


-0.078 


-0.192 


-0.443 



Table 8: Proton asymmetry coefficient in various scenarios. 





Oset-Saip| 


TRIUMFfSg 


UMass 


Y^T NM 77 (no corr.) 




1.6 


3.4 


it (corr.) 


1.5 


2.0 


0.5 
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1.2 


0.2 


^NM{p/n)n (no corr.) 




5.0 


4.6 


7i (corr.) 




5.0 


5.0 


n + K^; 7T, V , P ,U,K,K* 







4.0 


1.2 


r NM (PV/PC)ir (no corr.) 




0.4 


0.1 


7T (corr.) 




0.5 


0.1 


7T + A"[20|; n,7i,p,w,K,K*[U\ 




0.3 


0.8 



Table 9: Calculated properties of nonmesonic hypernuclear decay of j^C. 

3 — Hypernuclear Decay in Finite Nuclei 

Although the nuclear matter calculations are of great interest in identifying 
basic properties of the decay process, true confrontation with experiment re- 
quires calculations involving the finite nuclei on which the measurements are 
conducted. Of course, such calculations are considerably more challenging 
than their nuclear matter counterparts and require A shell model considera- 
tions as well as non-S-wave capture. Nevertheless a number of groups have 
taken up the challenge. For the case of the nonmesonic decay of ^ 2 C the re- 
sults are summarized in Table 9. In comparing with the experimental results 
given in Table 1, we see that the UMass calculation is certainly satisfactory, 
but the discrepancy between the UMass and TRIUMF work is disturbing 
and needs to be rectified before either is to be believed. 

A second nucleus on which there has been a good deal of work, both 
experimentally and theoretically, is |He, which is summarized in Table 10. 
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Oset-SalH 


TRIUMFpg 


TTBpll 


UMass 14 


^T NM TT(no corr.) 




1.0 


0.5 


1.6 


7r(corr.) 


1.15 


0.25 


0.144 


0.9 


7T + irpg ; 7T,r),p,u,K,K* 14] 




0.22 




0.5 


^NM{p/n)TT (no corr.) 




5.0 




15 


7i (corr.) 




4.8 




19 


tt + K 20 ; n,ri,p,u,K,K*\14\ 




5.4 




2.1 



Table 10: Calculated properties of the nonmesonic decay of |He. 

Here again what is important is not so much the agreement of disagreement 
with experiment but rather the discrepancies between the various calculations 
which need to be clarified before any significant confrontation between theory 
and experiment is possible. 

Before leaving this section, it is important to raise an additional is- 
sue which needs to be resolved before reliable theoretical calculations are 
possible — that of the AI = | rule.[]22| Certainly in any venue in which it 
has been tested — nonleptonic kaon decay — K — > 2n, 3n, hyperon decay— 
B — > B'n, AI = ~ components of the decay amplitude are found to be a 
factor of twenty or so larger than their AI = | counterparts. Thus it has 
been natural in theoretical analysis of nonmesonic hypernuclear decay to 
make this same assumption. (Indeed without it the already large number of 
unknown parameters in the weak vertices expands by a factor of two.) How- 
ever, recently Schumacher has raised a serious question about the correctness 
of this assumption, which if verified will have serious implications about the 



direction of future theoretical analyses. |23[ The point is that by use of very 
light hypernuclear systems one can isolate the isospin structure of the weak 
transition. Specifically, using a simple delta function interaction model of the 
hypernuclear weak decay process, as first written down by Block and Dalitz 
in 1963 one determines 

4 tt -n /„ /_A _ 2R nO 



yRe : 74 = T NM (n/p) 



3.Rpi + R p o 

A He : 7s = T NM {n/p) = — 

orCpi + rtpo 
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where here Rxj indicates the rate for N-stimulated hypernuclear decay from 
an initial configuration having spin j. One can then isolate the ratio R n o/R p o 
by taking the algebraic combination 



RnO 774 



R P o 1 + 74 - 775 



(20) 



and from the experimental values p5 



74 = 0.27 ±0.14, 75 = 0.93 ±0.55, 7 = 0.731^ (21) 
we determine 

R n0 0.20 

-0.12 



iL n 0.59T - 80 



(22) 

i p u.uy_ .47 

in possible conflict with the AI = | rule prediction — R n0 /R p0 = 2.[] If con- 
firmed by further theoretical and experimental analysis this would obviously 
have important ramifications for hypernuclear predictions. However, recent 
work at KEK has indicated that the correct value for 7 should be nearer to 
unity than to the value 0.73 used above in which case the ratio is consider- 
ably increased and there may be no longer any indication of AI = | rule 
violation. |26] 



4 — Conclusions 

We have given a brief overview of the field of weak hypernuclear physics. 
Because of limited experimental data and of the difficulty of doing reliable 
theory, the present situation is quite unsatisfactory. Although there is very 
rough qualitative agreement between theoretical expectations and experi- 
mental measurements, it is not clear whether discrepancies which do exist 
are due to experimental uncertainties, to theoretical insufficiencies, or both. 
On the theoretical side, what is needed are reliable calculations on finite hy- 
pernuclei (preferably by more than one group) which clearly indicate what 

1 Note that final state nn or np configurations which arise from initial 1 So states are of 
necessity 1=1. 
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signals should be sought in the data. The issue associated with the validity of 
the Al = ~ rule must be clarified. In addition there have been recent spec- 
ulations about the importance of two-nucleon stimulated decay[^7]] (which 
could account for as much as 15% of the decay amplitude according so some 
estimates) and of the importance of direct quark (i.e. non-meson-exchange) 
mechanisms, |28[] which deserve further study in order to eliminate the vexing 
double counting problems which arise when both direct quark and meson 
exchange components are included. On the experimental side, we require 
an extensive and reliable data base developed in a variety of nuclei in order 
to confirm or refute the predicted patterns. Clearly the strong program of 
hypernuclear physics at DA$NE will provide a major step in this direction. 
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